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The "14-1-11" phase compounds based on magnetic Mn ions and typified by 
Cai4MnBin and Bai4MnBin show unusual magnetic behavior, but the large number (104) 
of atoms in the primitive cell has precluded any previous full electronic structure study. 
Using an efficient, local orbital based method within the local spin density approximation to 
study the electronic structure, we find a gap between a bonding valence band complex and 
an antibonding conduction band continuum. The bonding bands lack one electron per for- 
mula unit of being filled, making them low carrier density p-type metals. The hole resides in 
the MnBi4 tetrahedral unit and partially compensates the high spin d 5 Mn moment, leaving 
a net spin near 4 fiB that is consistent with experiment. These manganites are composed of 
two disjoint but interpenetrating 'jungle gym' networks of spin | MnBi4 9 ~ units with ferro- 
magnetic interactions within the same network, and weaker couplings between the networks 
whose sign and magnitude is sensitive to materials parameters. Cai4MnBin is calculated to 
be ferromagnetic as observed, while for Ba^MnBin (which is antiferromagnetic) the ferro- 
and antiferromagnetic states are calculated to be essentially degenerate. The band structure 
of the ferromagnetic states is very close to half metallic. 
PACS: 71.20.-b,71.20.Be,75.30.-m,71.15.Mb 



I. INTRODUCTION 

Recently many new magnetic phenomena hav£ 
been discovered, such as colossal magnetoresistance 
(CMR), a new type of heavy fermion system 
(L1V2O4), and spin Peierls ground states in mag- 
netic insulators. Spin glass behavior has been found 
and studied intensively An magnetic insulators with- 
out structural disorderB. A related aspect of com- 
plexity in crystalline magnets arises when magnetic 
ions are distributed regularly but interionic dis- 
tances are large, ~1 nm. One example is the class 
of metallic rare earth (RE) silicides (RE)3Pd2oSi6, 
and the corresponding germanides, which present 
multiple ordering transitions with temperaturesQ in 
the range of a few K. Examples of dilute magnetic 
insulators are Naat^LiaF^u and the rare earth 
phosphomolybdatesQ (RE)P0 4 (Mo0 3 )i 2 x 30H 2 O, 
which order below 1 K, sometimes much below. A 
more well known example is the heavy fermion su- 
perconductor UBei3, which superconducts at 0>fl 
K and is near an antiferromagnetic instability!! 
These intriguing materials, with their complex ex- 
change interactions transmitted through intermedi- 
ate bonds, are examples of the continuing, perhaps 
even accelerating, growth and study of more com- 
plex compounds. The complexity may be structural 
in origin (many atomic sites of low symmetry in the 
unit cell), or it may be rooted in complex magnetic 



and electronic interactions (many different exchange 
couplings, perhaps competing and frustrating). 

In this paper we present a first-principles density- 
functional study of the electronic properties of two 
representative compounds of the so-called "14-1- 
11" phases, which are a clear example of materi- 
als showing simultaneous complex atomic and mag- 
netic structures. The 14-1-11 phasesEl typified by 
Cai4MnBin, where Ca may be substituted by Sr 
or Ba and Bi can be substituted by Sb or Asja are 
a rather difficult case of a magnetic system with a 
complex crystal structure. The structure will be 
described in detail below, but with its four for- 
mula units (104 atoms) and nine inequivalent sites, 
and magnetic ordering at 15-70 K indicating ex- 
change coupling of the order of 10 meV, this class 
presents a strenuous test for state of the art elec- 
tronic structure methods. Understanding the exis- 
tence of such a structure is itself an interesting topic 
in solid state chemistry. It seems to be consistent 
with the simplified-picture provided by the Zintl- 
Klemm-Bussmannt3 concept, which we will refer to 
as the Zintl concept. This picture represents a gen- 
eralization of the octet rule for binary semiconduc- 
tors and insulators, and invokes charge balance be- 
tween (nominally closed shell) structural units which 
themselves may be covalently or ionically bonded 
complexes as well as simple ions. However, both of 
the compounds that we address in detail in this pa- 



1 



per are metallidiinij in their magnetically ordered 
phases, whereas the Zintl arguments that have been 
applied to these compounds would be more appro- 
priate if they were semiconducting. In fact, the mag- 
netic ordering and the metallicity seem to be cor- 
related in these materials, which would be broadly 
consistent with a Rudennan-Kittel-Kasuya-Yosida 
(RKKY) type of modelEil for the origin of the mag- 
netic interactions. 

The two materials chosen for the present study 
are CauMnBin, which is ferromagnetic (FM) with 
Curie temperature Tc=55 K, and antiferromag- 
netic (AFM) Ba^MnBin, with Neel temperature 
Tat = 15 K. The observed magnetic moment of Mn 
and the ordering temperatures (several tens of K) 
give rise to several fundamental questions: 

(1) what is the charge state of Mn, and how does it 
relate to structural stability and conduction behav- 
ior? 

(2) given that the Mn ions are magnetic, why is the 
ordering temperature of the order of 50 K when the 
distance between Mn ions is at least 10 A? 

(3) what is the bonding path that provides the mag- 
netic coupling? 

(4) how is the magnetic order coupled to the carriers 
in these compounds, some members of which show 
colossal magnetoresistance near Tc? 

The only previous theoretical work on the elec- 
tronic structure of this class of compounds was done 
by Gallup, Fong, and Kauzlarich,o who considered 
a single-formula-unit simplification of semiconduct- 
ing CauGaAsn. They concluded that the bond- 
ing in this compound is consistent with the Zintl 
concept of valence counting using covalently bonded 
subunits (discussed below). The greatest difference 
between this compound and those we study in this 
paper is the substitution of the sp metal atom Ga 
in the site of inversion symmetry with the transi- 
tion metal atom Mn. Ga is trivalent, whereas the 
metallicity of the Mn-based compounds indicates a 
different valence for Mn. 

We conclude, in fact, that Mn is in a divalent 
state, with five 3d electrons localized and magnetic. 
This difference, trivalent Ga versus divalent Mn, 
leaves one unoccupied bonding orbital in the va- 
lence bands, giving metallic behavior as observed. 
The FM metal CauMnBin is close to a half metal- 
lic fillingEj of the bands. Our results suggest that 
adding one additional electron per formula unit to 
the itinerant valence bands should lead to a semi- 
conducting compound. Adding less than one more 
carrier could produce a half metallic FM situation. 
The most direct way of doing so would be the re- 
placement of a fraction of alkaline earth atoms with 
a trivalent atom, viz. Ca— >Y. A less likely possibil- 



ity would be Mn^Fe, if Fe would assume a trivalent, 
high spin d 5 configuration. The alternative is that Fe 
would assume a d 6 configuration, and therefore also 
be divalent, and the system would remain metallic. 

Current interpretation assigns aj-3j4n valence to 
the Mn atoms in these compoundai3'EJ. This dis- 
agrees with our findings as stated above, but it 
seems a very reasonable assumption considering that 
the measured magnetic moment is very close to 

4 hb /Mn for this class of materials. However, Mn 3+ 
also seems to imply a semiconducting character {viz. 
Cai4GaAsn), in contradiction to the experimental 
evidence. While the valence (or charge state) of an 
ion in a solid is a very useful concept, it does not nec- 
essarily represent an actual ionic charge; indeed, it is 
widely recognized that 'ionic charge' is an ill defined 
concept, particularly so for metals and narrow gap 
semiconductors. Nevertheless, the 'charge state' of- 
ten continues to be meaningful, and the charge state 
of the Mn atom in these materials has been deduced 
from measurements of the high temperature mag- 
netic suceptibility, which is indicative of moments of 
^4 fiB- This moment was assigned to Mn, arriv- 
ing at a d (trivalent) assignment. We find however 
that the density of states (DOS) contains a peak, 
well below the Fermi energy, originating from five 
bands per formula unit ( i.e. per Mn atom) which 
derive from Mn 3d states. This situation allows us 
to identify the Mn ion as having a d 5 configuration, 
consistent only with a divalent charge state. Nev- 
ertheless, the moment that we obtain is consistent 
with experiment, being closer to 4 /is than to the 

5 hb expected from a simple d 5 configuration. 
These two apparently contradictory observations 

are reconciled by the near-half-metallic band struc- 
ture: the holes (one per Mn atom) left in the valence 
band reside in the Bi4 tetrahedron that encapsulates 
the Mn ion. The effective ~4 \ib magnetic moment 
is still relatively well localized. As we will see be- 
low in more detail, both the experimental moment of 
~4 fxs and our observation of an almost half metallic 
band structure can be already anticipated from the 
electronic structure of the isolated MnBi^ " tetra- 
hedron, using the charge derived from applying the 
valence counting rules (n=9). 

The paper is organized as follows. Sections II 
and III describe the crystal structure and method 
of calculation, respectively. Magnetic energies and 
their interpretation in terms of exchange couplings 
are presented in Sec. IV, and Sec. V presents an 
analysis of the magnetization and its distribution 
along the unit cell. The density of states of crys- 
talline Cai4MnBin and Ba^MnBin are analyzed 
in Sec. VI, followed by a discussion of charge trans- 
fers and their relation with the formal valences in 
Sec. VII. Sec. VIII is devoted to examine the elec- 
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tronic states near Fermi level along with their re- 
lation with the anysotropic nature of the magnetic 
couplings. In Sec. IX the main results for the solid 
are reinterpreted as simple consequences of the elec- 
tronic structure of the isolated MnBi4 charged clus- 
ters, and a summary is presented in Sec. X. 



II. CRYSTAL STRUCTURE AND ITS 
IMPLICATIONS 

The structure of the alkaline earth metal pnic- 
tide compounds Ae^MPnn, which we will call the 
AeMPn structure, has been experimentally deter- 
mined, and descrihed in detail by Kurotomo, Kau- 
zlarich, and Webb.Eil In the present calculations we 
have used the coordinates given by these authors. 
Aci4MPnn compounds have a body center tetrag- 
onal (space group I4i/acd) unit cell with 4 for- 
mula units (104 atoms), a=17.002 A, c=22.422 A 
for CauMnBiii, and a=18.633 A, c=24.340 A for 
BauMnBiiiO The structure can be viewed as con- 
sisting of; interstitial alkaline earth atoms, isolated 
Bi (Bi3t3) atoms, distorted MnBi 4 (Bi2) tetrahe- 
dra, and Bi 3 (Bil-Bi4-Bil) linear units. The MnBi 4 
tetrahedra are translated by | along the c axis alter- 
nating with the Bi3 anions which are rotated by 90° 
with respect to each other, as shown in Fig. |l| The 
isolated Bi atoms are situated between the Bi3 and 
MnBi4 groups, along a screw axis which coincides 
with the c axis. All the Mn atoms in the unit cell 
are symmetry equivalent. The alkaline earth cations 
occupy four inequivalent sites, but the distinctions 
will not concern us. 

The MnBi4 tetrahedra are slightly flattened 
within the a-b plane, with the distortion increasing 
by ~1° when Ca ion is substituted by Ba (the in- 
equivalent angles are 118.0° and 105.4° for the Ca 
compound). All bond lengths in the system also 
increase with this substitution, for which a and c 
increase by 9-10%, with a corresponding volume 
change of nearly 30%. The Mn-Bi bond distance, 
for example, increases from 2.814 A (Ca) to 2.935 A 
(Ba). The Bi-Bi bond length in the Bi3 units also 
increases from 3.335 A (Ca) to 3.498 A (Ba)W, indi- 
cating that both this unit and the MnBi4 tetrahedra 
are substantially environment dependent. 

There are four cations located close to the Bi2 
atoms in the tetrahedron, with Bi-Ca distances in 
the range 3.2-3.3A, and 3.5-3.6 A for Bi-Ba sepa- 
rationscJ. The Bi4 central atom in the Bi3 unit 
is surrounded by four cations, while the terminal 
Bil atoms are coordinated with eight cations (Ca-[B|i 
bonds of 3.2-3.4 A, Ba-Bi distances of 3.5-3.75 A)B. 

As pointed out in Ref. [16], and will prove to be im- 
portant to understand the magnetic couplings, this 




FIG. 1. One view of the crystal structure of 
Cai4MnBiii, emphasizing the two interpenetrating sub- 
lattices of MnBI 4 tetrahedra, one shown as white and 
the other one darker. The Bi3 "sticks" are also shown; 
the Ca atoms and the isolated Bi atoms are not shown, 
for clarity. 



crystalline structure can also be regarded as two in- 
terpenetrating networks formed by MnBi4 tetrahe- 
dra. Each tetrahedron only belongs to one of these 
networks, and is linked to the four nearest MnBi4 
groups in the same network along the tetrahedral 
directions, i.e. through a Bi-Bi bond. This loose 
bond, with the Bi atoms 4.5-5 A apart, is medi- 
ated by three cations (atoms Al, Al', and A3 in 



Ref. 11), which are common nearest neighbors to 



both Bi (Bi2) atoms. Fig. || shows a schematic rep- 
resentation of this connectivity on the (100) plane. 

In contrast, the identification of the interaction 
pathway between MnBi4 tetrahedra belonging to dif- 
ferent networks is not so clear. Different paths can 
be envisioned, all involving several cations and at 
least one of the isolated Bi atoms (Bi3). 

For compounds such as Aei4GaPnn, formal va- 
lence arguments and accumulated experience indi- 
cate that the electronic structure can be rational- 
ized in terms of a model where each alkaline earth 
atom cedes two electrons (14A 2+ ), three of these 
valence electrons are collected by each isolated Pn 
atom (4Pn 3 ~), seven electrons can be transferred 
to the P113 units leaving unoccupied the antibond- 
ing a molecular stateE-3, and the remaining elecrons 
go to the tetrahedron centered on the Ga atom 
which becomes GaPn^. Although these formal 
charges are quite large and should not be taken lit- 
erally, such valence counting has proven to provide 
a very good description of the electronic structure of 
CaMGaAsnlii 

The metallic or semiconducting character of the 
Aei4MnBin materials therefore seems to depend in 
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TABLE I. Core radii (in a.u.) used for the generation 
of the pseudopotentials 





Mn 


Bi 


Ca 


Ba 


r s 


2.00 


2.20 


3.20 


3.50 


i P 


2.20 


2.90 


3.30 


4.00 


*d 


1.90 


2.90 


3.00 


3.50 



an essential way on the valence of the Mn atom. A 
trivalent Mn atom is isovalent with Ga and implies 
a semiconducting compound (which the AeuMnBin 
compounds are not), while a divalent Mn ion will im- 
ply a metallic compound (as observed) but appears 
to violate formal charge neutrality. One of the main 
objectives of this paper it to resolve this conundrum. 



III. METHOD OF CALCULATION 

The calculations-, ha.ve been performed with the 
program SIESTA,li3~EII an optimized code which al- 
lows standard density-functionalE3 calculations on 
systems with hundreds of atoms. This method 
has been successfully applied to the study of the 
electronic and_structural properties of mami-diffcr- 
ent materialscJ, including magnetic clustersEJ. The 
computational cost and memory requirements to 
build up the Hamiltonian matrix scale linearly with 
the size of the system, i.e. is of_pxder N (O(N)) 
where N is the number of atomsJlj'EEl while for the 
evaluation of the density and energy one can choose 
between standard methods, or make use of the re- 
cently developed 0(N) techniques .c3 In this work we 
have used a standard diagonalization of the Hamilto- 
nian because we are interested in the band structure 
and characteristics of specific states. 

The basis set is-||a_iinear combination of pseu- 
doatomic orbitals .OE§EZI In the ppsent-ealculations 
we have used a double-C polarizedEirE3~E3 (DZP) ba- 
sis set for all the atoms. A shell of d orbitals was also 
included for the Ba and Ca species. This amounts 
to 15 orbitals (2 s, 10 d, and a polarization (P) shell 
with 3 p orbitals) for Mn, 13 orbitals ( 2 s, 6 p, and 
5 P d) for Bi, and 10 orbitals ( 2 s, 5 d, and 3 P p) 
in the case of Ca and Ba, making 1192 orbitals in 
total. 

The core electrons are j-, re- 

placed by norm-conserving pseudopotentialstH gen- 
erated from the atomic configurations [Ar]3d 5 4s 2 for 
Mn, [Hg]6p 1/r5 6d - 25 5/ - 25 for BiB and, [Ar]4s 1 and 
[Xe^s 1 for Ca and Ba respectively. The core radii 
used in the generation of the pseudopotentials can 
be found in Table ffl. We apply the pseudopotentials 
using therj-ully separable formulation of Kleinman- 
By lander Ej 



The calculations have been carried out_in the lo- 
cal spin density approximation (LSDA)Lj, and a 
partial-core_correction for the non-linear exchange- 
correlatior£3 has been included for all the species. 
This correction is specially important to get reliable 
moments in magnetic atoms, but some care has to be 
taken in choosing the pseudocore radius. We found 
that a pseudocore radius of 0.70 a.u. leads to very 
good Jesuits in comparison to all-electron calcula- 
tions .L3 



IV. RESULTS FOR TOTAL ENERGIES AND 
MAGNETIC COUPLINGS 

In this section we report on the relative stabil- 
ity of FM and AFM alignments of the Mn moments 
and the magnetic couplings deduced from the cor- 
responding values of the total energy in both com- 
pounds. Since the calculations are already very large 
with the primitive crystallographic cell, we have not 
considered any magnetic alignments that would en- 
large the cell. Even so, with four magnetic Mn atoms 
in the cell there are three distinct types of AFM 
alignments, shown in Fig. ||in the (100) plane. The 
configurations are better explained in terms of two 
distinct ...-Mn-Bi-Bi-Mn-... 'chains' in this plane: 
(i) in-phase AFM chains (Fig. ||(a)), (ii) antialigned 
FM chains (Fig. 0(b)), and (iii) out-of-phase AFM 
chains (Fig. ||(c))! Similar chains run out of this 
plane and ultimately the various bonding 'chains' 
form two interpenetrating three dimensional net- 
works, as described in Section ||. Only four of the 
eight nearest Mn neighbors of a given Mn atom be- 
long to the same network as the central one, and only 
in the configuration (i) (in-phase AFM chains) the 
eight neighbors are antiferromagnetically aligned. 

From the figure it is clear that there are two dis- 
tinct nearest neighbor interactions, labelled Ji and 
.][ . The second neighbor interaction J 2 is also shown. 
We will obtain an estimate of these couplings assum- 
ing an effective spin coupling 

t <ij> 

where Sf = ±1 represents a normalized spin vari- 
able. The index t labels the types of couplings, and 
< ij > denotes spin pairs of that type. Note that 
in some treatments, our J might correspond to JS 2 
where S is the spin of the magnetic atom. 

For both compounds we have calculated the ener- 
gies for FM order and for the three AFM configura- 
tions decribed above. Of the AFM alignments, the 
antialignment of FM chains, as shown in Fig. ^(b), 
is always the one favored. Therefore, the minimum 
energy AFM configuration occurs when one of the 
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c-axis 



FIG. 2. Illustrations of the three inequivalent types 
of AFM order possible without enlarging the unit cell. 
A (100) plane is shown. The open circles with arrows 
indicate Mn with its spin direction, the filled dots indi- 
cate Bi sites, and the grey lobes indicate bonding pat- 
terns schematically. Some of the exchange constants are 
shown. 



TABLE II. Upper: Energies of the three antiferro- 
magnetic alignments pictured in Fig. ^ relative to the 
ferromagnetic alignment. The type of AFM order in 
Bai4MnBin (type (a), (b), (c), or other, is not estab- 
lished experimentally). Numerical uncertainties are ex- 
pected to be ~ ± 5 meV. Lower: Exchange couplings 
for a short range Ising model determined from the total 
energy differences. The couplings correspond to those 
pictured in Fig. ^. J2+3 are the value of the couplings 
if second and third neighbors would enter in the energy 
expression via a single exchange constant (see the text). 



Energy (meV) 


Cai 4 MnBin 


Bai 4 MnBin 


AFM(a) - FM 


59 


12 


AFM(b) - FM 


25 


4 


AFM(c) - FM 


52 


13 


Experiment 


FM 


AFM 


Exchange Constants (K) 


Cai 4 MnBin 


Bai 4 MnBin 


Ji 


63 ± 5 


15 ± 5 


A 


23 ± 5 


2 ± 5 


■h 


27 ±10 


7 ±10 




9 ± 3 


2 ± 3 


Tc, mean field 


130±20 K 


28±20 K 


Tc, experiment 


T c =55 K 


Tjv=15 K 



Mn networks is entirely spin up, the other entirely 
spin down. 

Since the energy differences are quite small (as ex- 
pected), it was necessary to study the convergence 
of the energies-.with respect to both the number of 
points (cutofO) in the real-space grid, and the It- 
point mesh used to sample the Brillouin zone. Cut- 
offs up to 216 Ry for CauMnBiu and 147 Ry for 
Bai4MnBin were used, with up to 12 k-points in 
the irreducible Brillouin zone. The results shown in 
the following have been obtained with 3 inequiva- 
lent k-points, and a cutoff of 150 Ry and 147 Ry for 
the Ca and Ba compounds, respectively. With these 
parameters the total energies are converged, for a 
given basis set, to 10 meV/Mn, and 5 meV/Mn for 
the energy differences between different configura- 
tions, which show a faster convergence. 

The calculated total energy differences are pre- 
sented in Table || For CauMnBiu the FM state 
is more stable by 25±5 meV per Mn atom than the 
lowest energy AFM phase, in accord with the ob- 
servation of ferromagnetism in this compound. In 
a simple nearest neighbor Ising or Heisenberg mod- 
elling of the spin coupling this energy would corre- 
spond to a FM exchange coupling of J^3 meV. For 
Bai4MnBin the energy difference is much smaller, 
with the FM phase 4±5 meV/Mn more stable, i.e. 
degenerate to within our accuracy. This compound 
is observed to be AFM, with Tn—15 K, with the low 
ordering temperature reflecting smaller magnetic in- 
teractions than in Cai4MnBin. 
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Tabic H also presents the exchange couplings of 
the effective spin Hamiltonian (hi). Each Mn atom 
has eight Mn first neighbors, four on the (100) and 
four on the (010) planes, at distances of 10.18 A for 
the Ca compound (10.84 A for Ba^MnBin). Al- 
ready from the energies in Table |J, it is evident that 
two different first neighbor couplings need to be in- 
cluded. In fact, if just one exchange constant were 
to be used in a model with only first neigbor inter- 
actions, the energies of configurations (b) and (c) in 
Fig. |^ should be identical. Physically, the need of 
two constants is easily understood as a consequence 
of the presence of two Mn networks: stronger in- 
teractions can be expected within the same network 
(Ji) than between atoms in different networks (J^). 
This is clearly confirmed in Table [ill where Ji is 
shown to be much lager than for both compounds. 

There are two second neighbors located at 11.21 A 
(12.17 A) along the c axis. The four third neigh- 
bors can be found 12.02 A (13.18 A) away, along 
the <100> and equivalent directions. In the crys- 
tallographic cell these last atoms are equivalent to 
the second neighbors yet located at somewhat larger 
distances due to the tetragonal distortion of the cell. 
Further neighbors are at more than 17 A for both 
compounds and their interaction constants are ex- 
pected to be much smaller. With the calculated en- 
ergy differences we cannot make independent esti- 
mations of the second (J2) and third (J3) nearest 
neighbor couplings. At this point, we can choose to 
restrict the interactions up to second neighbors and 
obtain a value for J2. Another possibility is to con- 
sider second and third neighbors as entering in (4.1) 
at the same footing, i.e. through an averaged cou- 
pling constant (J2+3=(2J2+4J 3 )/6). The values for 
these exchange constants are also listed in Table |lj]. 



V. MAGNITUDE AND DISTRIBUTION OF 
THE MAGNETIC MOMENT 



In this section we study the distribution of the 
magnetic moment in the unit celLwith the help of 
the Mulliken population analysisa. It must be rec- 
ognized that the populations obtained using Mul- 
liken analysis suffer, as those obtained via all simi- 
lar techniques, from an inherent arbitrariness: both 
the charge density and the total charge are observ- 
ables of the system, but their partition into differ- 
ent atomic-like contributions cannot be uniquely de- 
fined. Accordingly, Mulliken population trends and 
differences are more meaningful that their absolute 
values, and their use is fairly common. We obtain 
in this way information about the effective valence 
state and magnetic moment of the Mn atoms, the 
polarization of the surrounding atoms, and the de- 



TABLE III. Population and moments of the Mn 
d shell, and magnetic moments of Mn atoms and 
MnBi4 tetrahedra as obtained from a Mulliken popula- 
tion analysis for the FM Cai4MnBin, and Bai4MnBin. 
HAexiMnBin stands for the net magnetic moment per 
formula unit, while A^// * s the effective moment per Mn 
atom as obtained from the measurements of the high 
temperature magnetic suceptibility. 





Cai 4 MnBin 


Bai 4 MnBin 




4.67 


4.74 


05 


0.51 


0.42 


/id 


4.16 


4.32 




4.45 


4.62 


HMnBi 4 


4.22 


4.35 


/-lAei^MnBin 


4.25 


4.40 


exp 
Veff 


4.8 


4.8 



gree of localization of the total magnetic moment. 

The magnetic moments on the Mn atoms for the 
FM compounds, as obtained from the difference 
between the Mulliken populations of the majority 



and minority spin densities, are shown in Table III. 
They are found to be nearly independent of the type 
of magnetic order, and differ by 4%: 4.45 hb for 
CauMnBiu, 4.62 fi B for BauMnBiu. The part of 
the moment attributable to the 3d shell is 4.16 [is 
and 4.32 /xg, respectively, with about 0.1 [Ib and 
0.2 us induced in the 4s and Ap states in each com- 
pound. The spin up population (Table III) of the 



Mn 3d orbitals approaches 5 electrons, indicative of 
a d 5 configuration. Several factors can contribute to 
a reduction of the 3d population, being the most im- 
portant the hybridization with the 6p states of the 
neighboring Bi atoms. The identification of the con- 
figuration of Mn as 3d 5 will be also confirmed from 
a detailed a nalys is of the electronic band structure 
(see Section VI B ), providing further support to our 
picture of a divalent Mn atom. 

The total moment per formula unit of the FM 
state is calculated to be (see Table III) 4.25 hb for 
CauMnBiu and 4.40 /j,b for Ba^MnBin, which 
becomes 17 \xb and 17.6 /j,b per cell, respectively. 
With all the majority d states occupied, each Mn 
atom should contribute with a magnetic moment of 
+5 fj,B (+20 [lb in the unit cell). However, analysis 
of the contributions to the d bands indicates that 
this moment, although still mostly localized in Mn, 
is somewhat spread out to the neighboring Bi atoms 
(Bi2) due to the considerable hybridization between 
the Mn 3d states and the 6p states of Bi. Therefore, 
it is more safely pictured as associated with the five 
atoms in the MnBi 4 tetrahedra, rather than to the 
Mn atoms alone. 

A reverse polarization of 2.4-3 /is (0.6-0.75 
/is/MnBi4 unit) must be present to recover the cal- 
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culated moment. From the data in Table [II it is 



clear that this back polarization is also concentrated 
in the MnBi4 tetrahedra, whose magnetic moment 
almost recover the total value per formula unit. Part 
is due to the occupation, through hybridization, of 
the minority spin Mn d orbitals. The remainder is 
mostly divided among the Bi2 atoms, which have a 
magnetic moment of ~ —0.06 lib per Bi atom, the 
largest in the cell after Mn. This moment is quite 
modest, but it is in fact built up from two much 
larger contributions that tend to cancel each other: 
~ +0.20 [ib coming from their participation in the 
majority Mn 3d bands (as previously commented), 
and a ~ —0.26 lib component coming from the po- 
larized Bi 6p bands reponsible for the reverse po- 
larization. The amount of spin polarization spread 
around the rest of the cell is negligible (~ 0.001 lib 
per Bi or Ae atom on average). Only the isolated 
Bi atoms (Bi3), and the alkaline earth atoms (Al 
and A3) connecting the MnBi4 tetrahedra present a 
small polarization of ~ —0.03 lib and ~ +0.02 lib, 
respectively. 

An analysis of the Mulliken populations of all 
the studied AFM orders, for both Ca and Ba com- 
pounds, indicates a very similar picture: the calcu- 
lated moment per formula unit has to be assigned 
to the MnBi4 tetrahedra, being neither fully local- 
ized in the Mn ion nor spread over the whole unit 
cell. This moment, of ^Allb, is built up from two 
distinct contributions, a ^5lib moment coming from 
the fully occupied majority spin 3d shell of Mn, and 
a ~ — 1/j,b reverse polarization coming from the 6p 
states of the four Bi atoms in the MnBi4 complex. 
The MnBi4 tetrahedra are therefore the magnetic 
units in these compounds. 

Experimentally, an effective moment (/x e //) of 

4.8 lib has been obtained for both compounds from 
the high temperature suceptibilitytl In fact, for 
all the related compounds Ae14MnPn.11, where Bi is 
substituted by Sb or As, and the cation sites are oc- 
cupied by Ca, Sr or Ba atoma.,|-£w / is always found 
within the rangepA.Sg-5.4 nJMR For YbuMnBin 
and Ybi4MnSbnE3 , EJ Li e ff is also measured to be 

4.9 hb- This has been interpreted, using the stan- 
dard form for the Curie constant, /x^j=g 2 J(J+l)/x^ 
where J is the total effective angular moment and g 
is the Lande g-factorO, as a signature of the pres- 
ence of localized ^4 lib moments in all these com- 
pounds. This seems in good agreement with our cal- 
culated magnetic moment per formula unit (i.e. per 
M11B14 tetrahedron). However, all previous works 
have assigned these magnetic moments solely to the 
Mn atoms, leading to the conclusion that Mn is in 
a 3+ valence state with four unpaired d electrons, 
which is in contradiction by our findings. 

A better quantity to compare with our calculation 



of the magnetic moment is the experimental satura- 
tion moment, which is the value of the ordered mo- 
ment parallel to the applied field <M 2 >. TliCjiaitial 
work on hysteresis loops for powder samples, Ertil ob- 
served saturation moments of 2.5-3 /xs/Mn, which 
are ~30% smaller than both the estimations based 
on the measured Curie constants, and our calculated 
total moments. The polarization of the itinerant va- 
lence electrons was invoked in those works to explain 
this discrepancy between the measured saturation 
moments and the fitted Li e ff- Ho wever , more recent 
measurements using single crystalsHlill indicate that 
the saturation moment is close to 4 /xg/Mn, in bet- 
ter agreement with our calculations. 



VI. DECOMPOSITION OF THE 
ELECTRONIC SPECTRUM 

Besides the character of the band structure near 
the Fermi energy, which will be discussed in detail in 
a separate section, the general features of the elec- 
tronic structure and ordering of the levels are the 
same for both FM and AFM orders, and very sim- 
ilar for Ca and Ba compounds. Therefore we will 
concentrate here in the FM order. 



A. Total Density of States, General Structure 

With five valence electrons for Bi, two for Ca, and 
seven for Mn, and four formula units per cell, the 
occupied valence bands of these compounds must 
accommodate 360 electrons, or roughly 180 occupied 
bands of each spin. 

The total density of states (DOS) of FM 
Cai4MnBin and Ba^MnBin are shown in 
Fig. U (a) and (b), respectively, which also reveal 
some general features. These plots have been ob- 
tained using three inequivalent k-points and a gaus- 
sian broadening of 0.1 eV. Tests indicate that using 
more complete k-sampling (up to 12 inequivalent k- 
points) produces almost identical results. The oc- 
cupied valence bands are 4-4.5 eV wide and nearly 
filled. A 'gap' of ~ 1 eV separates the valence 
complex from a continuous spectrum of unoccupied 
states coming mainly from the Ca-Bi hybridization. 
Within the gap, 0.5 eV above the top of the va- 
lence band, we find four narrow (0.2 eV) bands. In 
the minority spin DOS, near the top of the gap for 
Cai4MnBin, and confused with the onset of the con- 
duction band in BauMnBin, we can also find a peak 
that contains the unoccupied Mn 3d bands. 

The band complexes below the gap consist of 182 
bands of each spin for AFM states, for which the 
spin directions are equivalent. For FM order there 
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are 192 majority and 172 minority bands, as the ma- 
jority bands contain twenty more bands (for the five 
occupied 3d electrons of the four Mn atoms) than 
the minority bands. The total number of bands 
in the valence complex is consequently the same 
for FM and AFM states, four electrons lacking for 
these compounds to become semiconductors. The 
172 bands not associated with the d states of Mn 
have mainly Bi character. 44 of them, located ap- 
proximately 9 eV below the Fermi energy, are due 
the Bi 6s states. The remaining 128, which form 
the uppermost valence band complex, can be clas- 
sified as Bi 6p bands, exhibiting some degree of hy- 
bridization with Mn and Ca/Ba states. From the 
44 Bi atoms present in the unit cell we could expect 
to have 132 Bi 6p bands, however. The four Bi 6p 
bands left out from the valence complex form the, 
already mentioned, narrow peak in the middle of the 
gap. These bands correspond to very localized states 
associated with the Bi3 chains (see below). 

The DOS at the Fermi energy (Ep) is higher for 
the majority spin, (Fig. ||) and both compounds 
are very close to half metallic, with the minority 
bands almost fully occupied (as discussed later). 
From the specific heat Siemens et alS3 obtained 
the Fermi level DOS N(E F ) = 1.7±0.7 states/eV- 
atom for Ba^MnBin, and 1.6±0.1 states/eV-atom 
for the similar compound Sr^MnBin. These are 
quite large values, comparable to the DOS of d 
band metals. However, Siemens et al. pointed out 



that contributions to the specific heat coming from 
the nuclear hyperfine splitting, not taken into ac- 
count in their analysis, could obscure their results, 
leading to large apparent values of the electronic 
DOS. In fact, our estimations of the DOS at Ep 
are much smaller: N(E F ) = 0.3 states/eV-atom for 
both Ca and Ba compounds with FM order. For 
AFM order, N(E F ) = 0.40-0.46 states/eV-atom de- 
pending on the specific AFM alignment. These val- 
ues can be compared with the average DOS in the 
valence band complex; leaving aside the Mn bands, 
we have 256 bands for both spin orientations dis- 
tributed over a range of ~4 eV, leading to an average 
<N(E)>=0.6 states/eV-atom. 



B. Mn Character 

Because of the magnetic moment on the Mn atom, 
we begin our discussion of the electronic bonding 
with the Mn d states. In all the cases we have consid- 
ered, the simple characterization is that the majority 
d states are filled and the minority states are empty, 
with an exchange splitting between them of ~4 eV. 
Accordingly, the Mn atom can be described as hav- 
ing a id 5 configuration or, alternatively, a 2+ va- 
lence. The MnBi4 tetrahedron is compressed along 
the a — b plane and only the d xz and d yz states are 
degenerate, so the crystal field results in four dis- 
tinct d levels. The crystal field splitting is negligi- 
ble, however. In Ca^MnBin (BawMnBiii), shown 
in Fig. |], the majority d states lie in a single peak 
centered 2.7 eV below Ep with a width of roughly 1 
eV (respectively 0.4 eV). The d states do hybridize 
with all states above this peak (with mainly Bi 6p 
character) to the gap just above Ep. The minor- 
ity states are unoccupied and are concentrated in 
a peak 1.3 eV above Ep. This unoccupied peak is 
narrow in Cai4MnBin (0.6 eV) and lies within the 
gap (see also Fig. ||). In Ba^MnBin it is much 
broader (at least 1 eV) meeting the bottom of the 
conduction band, and exhibiting a higher degree of 
hybridization with the conduction states. 



C. Bi Character 

The spin-averaged DOS projected into the differ- 
ent types of Bi atoms in the unit cell can be found in 
Fig. g. The Bi 6s bands lie 9-10 eV below the Fermi 
level, and have negligible hybridization with states 
of different character. The Bi 6p states make up a 
large component of the valence bands, exhibiting hy- 
bridization with both Mn and alkaline earth states, 
and hence are crucial in the bonding and in the mag- 
netic coupling. The valence bands are 4 eV wide in 
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Bai4MnBin and slightly (5%) wider in the smaller 
volume Cai4MnBin compound. There is a 'gap' of 
~1 eV between the valence and conduction bands, 
except that a very narrow set of unoccupied bands 
lies within this gap. Although we have neglected in 
Fig. |s| the small effects due to the induced exchange 
splitting, it should be kept in mind that the polar- 
ization of the Bi 6p states near the is crucial to 
understand the magnitude and distribution of the 
magnetic moment. 

The complications are still considerable, since 
there are four distinct Bi sites. The Bi2 sites corre- 
spond to the distorted MnBi4 tetrahedron, the B13 
linear chain is formed by a central Bi4 atom bonded 
to two Bil atoms and there are also four isolated 
Bi3 atoms. The first thing that is evident from the 
Bi spectra in Fig. || is that there is not a great deal 
of difference in the densities of states of the four Bi 
sites, indicating that their charges are not very dif- 
ferent. The Bil-Bi4-Bil linear unit has only a small 
(but nonzero) DOS at Ep, whereas both the Bi2 
bonded to the Mn and the 'isolated' Bi3 site have 
much larger fractions of the Fermi level DOS. In the 
case of the Bi2 atoms this DOS is almost entirely 
due to the majority spin contribution. 

The unoccupied 'gap states', one state and there- 
fore one band for each Bi3 unit (four in the unit 
cell), are associated with the Bil and Bi4 atoms 
only, with no appreciable spin polarization. Since 
they lie in the gap and have a dispersion of only 0.2 




Energy (eV) 

FIG. 5. Densities of states projected into the s (long 
dashed lines at lower energies), p (solid line), and d 
(dashed lines at higher energies) symmetry orbitals of 
the Bil, Bi2, Bi3, and Bi4 sites in ferromagnetically 
aligned Cai4MnBin. The differences occur primarily in 
the gap region 0-2 eV, and are discussed in the text. 
Spin-averaged spectra have been plotted; only the po- 
larization on the Bi2 site (~ — 0.06/is) is appreciable. 



eV, these states are quite localized. Fig. |6| shows 
contour plots of their density in a plane containing 
the Bil-Bi4-Bil sites. Their shape, with nodes in 
the bond region between Bil and Bi4 atoms, and 
the fact that (see Fig. ||) they have a larger weight 
in the central Bi4 atom, reveal them as antibonding 
er molecular orbitals. The presence of these a* states 
and the negligible population of the 5c? states of the 
Bi4 atom, which rules out the possibility of a dsp 3 
hybridization, confirm the description of the Bi3 7 ~ 
anion as a hypervalent "three center, four electron" 
bonded structure. This type of bonding was already 
predicted for AS3 in CauGaAsn by Gallup et a/.EJ, 
who applyed first-principles calculations to a simpli- 
fied model (30-atoms unit cell) of this compound. 

The conduction bands contain the rest of the Bi 
6p spectral weight along with Bi 5d, Ca or Ba s, and 
Mn 3s, 3p and minority 3d contributions. 



D. Ca/Ba Character 

The formal valence picture suggests considering 
the Ca and Ba atoms as dipositive ions. This picture 
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FIG. 6. Contour plot in the (001) plane of the density 
of the antibonding a* states localized in the Bi3 chains 
of FM Cai4MnBin . These four bands are located within 
the gap, around 1 eV above the Fermi energy. The con- 
tours begin at 3xl0 -4 e/Bohr 3 and increase by steps 
of 6xl0 -4 e/Bohr 3 (respectively, 7 and 14 electrons per 
unit cell). The positions of the Mn and Bi atoms in the 
plane are indicated schematically. 

should not be taken too literally, as we find some s, 
p, and d alkaline earth contributions to the DOS 
in the valence band region (figure not shown). The 
Ca/Ba d DOS peaks lie 4 eV above Ep, whereas the 
s and p contributions are spread over a large region 
from ~-3 eV through and above the d bands, with 
rather little structure in their distribution. 



VII. CHARGE TRANSFER AND FORMAL 
VALENCE 

Now we address the question of whether the for- 
mal valence that the Zintl picture associates to the 
different groups of atoms provides a reasonable de- 
scription of our results. The formal valences are: 
Ca 2+ /Ba 2+ , isolated Bi 3 ~, Big" linear chains, and 
MnBi^ - tetrahedra. These latter two are very 
large charges and, as can be expected, the Mul- 
liken populations are much smaller, although still 
indicative of large charge transfers. For example for 
Cai 4 MnBin we obtain: Ca 8+ , MnBi 4 3 7 -, Big 2 8 ", 
Bi 12- . The main reason for this difference between 
formal and Mulliken charges can be traced back to 
the fact that, while there is not a single state with 
pure or main Ca/Ba character, most of the states in 
the valence band exhibit some small hybridization 
with the orbitals of the alkaline earth atoms, spe- 
cially those at the top of the valence band with the 
d states of Ca/Ba. In fact, the s population of the 
alkaline earth ions (the only orbital that would be 
included in a simplified description of these atoms) 
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FIG. 7. (a) Contour plot of the deformation den- 
sity of ferromagnetically aligned Bai4MnBin in a (110) 
plane, illustrating the Bi3 unit bonding (top panel). The 
bottom panel (b) shows the (100) plane through the 
MnBi4 tetrahedron, indicating the charge 'deformation' 
is rather spherical around the atoms and is therefore 
primarily charge transfer. Only the region of increased 
charge is shown, with contours beginning at zero and in- 
creasing by 0.03 e/A 3 . The position of the atoms in the 
planes is indicated by their chemical symbol. The shaded 
areas correspond to regions with defect of charge. 



FIG. 8. Greyscale plot in a (100) plane of the 
(pseudo-)charge density associated to the unoccupied 
states below the gap (holes) in AFM Ca^MnBin . Panel 
(a) shows the spin up holes, and (b) the spin down holes. 
The arrows indicate schematically the positions of the 
Mn atoms and the orientations of their atomic magnetic 
moments. The positions of one of the Bi2 and Bi4 atoms 
are also indicated. We have used a logarithmic scale 
which saturates to black at 10 -3 e/Bohr 3 and to white 
at 10 -5 ' 3 e/Bohr 3 (respectively, 22 and 0.1 electrons per 
unit cell) 
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is only 0.4 electrons for both Ca and Ba compounds, 
picturing them as Ca(Ba) 16+ . Futhermore, Ca/Ba 
components are generally very small for each indi- 
vidual eigenstate of the solid, and formally they can 
be safely considered as dipositive. 

Fig. |t] shows the contour plots of the deformation 
density (difference between the self-consistent charge 
density and that obtained from the sum of spherical 
neutral atom contributions) in different planes for 
FM Bai4MnBin. Two things can observed: i) there 
is a great deal of charge transfer to the Bi atoms 
from the Ca/Ba atoms, which are always surrounded 
by a region showing a depletion of charge; ii) the 
deformation density is quite featureless and spherical 
around the Bi atoms, indicating that indeed the ionic 
charge transfer is the main mechanism involved. 

We now consider the polyatomic anions. The sit- 
uation of Bi3 is specially clear. Two antibonding 
molecular states (one per spin) remain unocuppied 
and the valence states are almost completely occu- 
pied up to the gap, with a very small contribution 
to the Fermi level DOS (see Figs. | and |). All these 
data together identify the Bi3 unit as having a 7- for- 
mal valence. The weight of the gap peaks in Fig. || 
suggest that less charge is transfered to the central 
Bi4 than to the Bil sites, which is confirmed by the 
Mulliken analysis (around 0.3 |e| less). This com- 
pletely agrees with the "three center, four eLeetrons" 
model of bonding proposed for the Bi3 unitllJ. 

We now discuss the cases of the MnBi4 units and 
the isolated Bi3 ions. On one hand, a 3- formal 
charge would correspond to closed shell Bi3 anions, 
with a negligible contribution to the DOS at Ep-. 
On the other hand, considering the 2+ effective va- 
lence (d 5 ) of the Mn atom, the MnBi4 group would 
be formally able to accept up to ten electrons. This 
implies that a MnBi^ - tetrahedron would be lack- 
ing one electron to be closed shell. This agrees 
with our analysis of the electronic structure (Sec- 
tion VIA), which indicates that these compounds 
lack four electron per unit cell (i.e. one electron 
per MnBi4 group) to become semiconductors. Fur- 
thermore, taking into account the calculated reverse 
polarization of the Bi2 sites, the hole associated to 
each MnBi^ - tetrahedron should be expected to be 
parallely aligned with the atomic moment of the cor- 
responding Mn. This is confirmed by the data pre- 
sented in Fig. ||, where a plot of the density associ- 
ated to the holes for AFM CauMnBin (alignment 
type (a) of those shown in Fig. ^J) in the (100) plane 
is shown. This plane contains Mn atoms, some of 
the Bi2 atoms bonded to them, and Bi4 sites. We 
can verify how the holes, while not exactly confined 
to, are fairly localized in the tetrahedra: the spin up 
holes in those tetrahedra where the atomic moment 
of Mn points in the up direction, and conversely for 
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spin down holes. An obvious consequence of this po- 
larization of the valence band holes is the (close to) 
half metallic character of the band structure of the 
FM compounds, which will be examined in detail in 
the next section. Hence we can conclude that the in- 
dentification of a valence 9- for the MnBi4 unit and a 
valence 3- for the Bi3 ions are meaningful, and from 
the point of view of a general picture our results are 
consistent with counting of formal charges. 



VIII. ELECTRONIC STATES NEAR THE 
FERMI LEVEL 

A. Character of the states 

For simplicity we will first discuss the case of 
the AFM order. Figj9| shows the band structure 
for AFM Cai4MnBin (alignment type (b) of those 
shown in Fig. [2]). There are two bands unfilled be- 
low the gap, and hence there are four electrons less 
(i.e. one per formula unit) than required to make 
this compound a semiconductor. The figure shows 
with detail the region of the 'gap' between bonding 
and extended conduction states. For each spin ori- 
entation this gap contains two groups of bands cor- 
responding to very localized states: one very narrow 
group of four bands coming from the antibonding 
states of the Bi3 chains in the energy range 0.85- 
1.05 eV, and ten unoccupied 3c? Mn bands in the 
range 1.2-1.6 eV, the last ones almost degenerate 
with the begining of the continuum of conduction 
states. The energy position of these groups of bands 
is very similar for AFM and FM states, but while 
the Bi3 bands do not show any spin polarization, of 
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FIG. 10. Band plots in the vicinity of the Fermi level 
along primary symmetry directions for ferromagnetically 
aligned Bai4MnBin and Cai4MnBin . The symmetry 
points correspond to the tetragonal BCC unit cell. The 
energies are referred to the Fermi level. 

course for FM order the twenty unoccupied Mn 3d 
bands only occur in the minority band structure. 



B. Half metallicity of FM compounds 

The bands for both FM Ca^MnBin and 
Bai4MnBin near Ej?, along several symmetry direc- 
tions, are shown in Fig. [lo| Except for details (which 
affect Fermi surface shape and therefore might be- 
come important) the bands are similar and we dis- 
cuss their common characteristics. The band struc- 
ture is very close to half metallicity. The minority 
bands are occupied except for a T-centered pocket 
that contains only a fraction of a hole (around 50%). 
The rest of the four holes correspond to the majority 
bands, which have several partially occupied bands: 
five in Cai4MnBin, four in Bai4MnBin. 

The almost half metallic band structure is 
especially interesting in light of the reported 
colossal magnetoresistance in, the related com- 
pound Euia.gyGdo.oaMnSbny The CMR system 
Lai_ x D x Mn03, with D = Ca, Sr, or Ba, and x ~ i, 
are believed to be half metallic or close to itEJ. 
Eui3.97Gdo.o3MnSbn represents a case in which a 
compound, Eu^MnSbn, that is isovalent with the 
ones studied here is doped with additional carriers 
(Eu 2+ — > Gd 3+ ), which will drive it toward half 
metallicity. 



FIG. 11. Greyscale plot in a (100) plane of the 
(pseudo-)charge density near Ef in FM Cai4MnBin . 
Top panel (a) for occupied states, and bottom panel 
(b) for holes. The positions of one of the Bi2, Bi4 
and Mn atoms in the plane are schematically indicated. 
Note the bonding chains coupling Mn atoms (darkest 
spots) formed by the two intervening Bi2 atoms. We 
have used a logarithmic scale which saturates to black 
at 10" 2 ' 8 e/Bohr 3 and to white at 10~ 5 ' 3 e/Bohr 3 (re- 
spectively, 35 and 0.1 electrons per unit cell) 

C. Magnetic Couplings 

An RKKY model for the magnetic interaction is 
the simplest assumption to make in a dilute mag- 
netic metal, and in the absence of other informa- 
tion, that is what has been used so far in the 
interpretation, of magnetic coupling in these 14-1- 
11 magnets.E3 Our identification of near neighbor 
couplings (Table ||) invalidates a spherical RKKY 
model, where the coupling strength depends only on 
distance. This conclusion is supported by Fig. [ll], 
which shows a grayscale plot of the (pseudo-)charge 
density from states near E^ in Cai4MnBin in a 
(100) plane containing -Mn-Bi2-Bi2-Mn- bonding 
chains. In this figure, panel (a) contains the den- 
sity from states in the range from E^-0.5 eV to Ep. 
Panel (b) contains the charge density from the unoc- 
cupied states between E^ and the gap (holes). Both 
panels clearly show that the density associated to 
the states near Ep concentrates along the -Mn-Bi2- 
Bi2-Mn- chains, making them the most likely con- 
duit for the exchange coupling between Mn atoms. 
The occupied states, panel (a), shows direct charge 
peaks along this bonding chain, although the states 
may be predominantly antibonding. Panel (b) also 
indicates density along this chain, and is apparently 
more antibonding, since the charge maxima centered 
on the Bi2 sites do not point toward the neigh- 
boring Mn atom. However, this has to be taken 
with some caution since, as already discussed in Sec- 
tion [n], other atoms may intervene in the bonding 
between neighboring MnBi4 units. In particular, 
three Ca/Ba common nearest neighbors of those Bi2 
atoms in neighboring connected tetrahedra, which 
are not contained in the (100) plane shown. In fact, 
the charge maxima in Fig. [ll] (b) seem to point from 
the Bi2 atoms to some of those alkaline earth ions 
located ~ ±0.7 A over the plane. 

The general picture of the magnetic couplings that 
arises from our calculations is the following: Mn 
atoms are divided between two different networks 
(defined by the -Mn-Bi2-Bi2-Mn- bonding chains); 
the nearest neighbors coupling within the same net- 
work is the strongest and always FM; the couplings 
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FIG. 12. Energies of the molecular orbitals of the 

MnBi^ - and MnBiJj" tetrahedral clusters. The zero of 
energy is located at the highest occupied molecular or- 
bital. The cluster geometry was fixed to that of the 
MnBi4 tetrahedron in crystalline Ba^MnBin. See text 
for notation. 



among Mn atoms belonging to different networks are 
much smaller and may be AFM, explaining the ex- 
perimental AFM ground state of Ba^MnBin. In 
particular, the presence of two internally FM cou- 
pled networks is consistent with the observation 
that, although it suffers an AFM transition, the 
Curie temperature obtained for Bai4MnBin from 
the fitting of theJaigh temperature magnetic sucepti- 
biliy is positive.Eil Therefore, although we have failed 
to predict the antiferromagnetism of Ba^MnBin, 
we expect that the magnetic alignment in this com- 
pound must be quite similar to the most stable of the 
AFM orders discussed in this paper (Fig. |^ (b)), i.e. 
one network entirely spin up, and the other entirely 
spin down. 



IX. DISCUSSION: STUDY OF CHARGED 
MNBI4 CLUSTERS 



The results of the previous sections point out 
that, rather than the Mn ions, the MnBi4 tetrahedra 
should be considered as the magnetic units in these 
compounds. The amount of charge transferred to 
the MnBi4, in combination with the effective valence 
of the Mn atom, also plays a crucial role to under- 
stand the magnetic moment, the conduction prop- 
erties, and the interplay between them. For this 
reason we devote a separated section to the study 
electronic structure of free standing charged MnBi 4 
tetrahedral clusters and how these results translate 
into implications for the solid. 



Fig. |U shows some of the energy levels of clus- 
ters with a total net charge of -9 |e| (nominal for 
the MnBi4 group in the Aei4MnBin compounds), 
and a smaller charge of -6 |e|. The geometry of the 
cluster was kept fixed to that of a MnBi4 tetrahe- 
dron in crystalline Ba^MnBin. Due to the use of 
periodic boundary conditions, it was necessary to as- 
sume a uniform neutralizing background in the cal- 
culations. The horizontal scale represents the pro- 
jection of each molecular orbital (MO) into the dif- 
ferent atomic species, solid and thinner dashed lines 
correspond to Mn and B14 contributions, respec- 
tively (which sum to unity). Different degeneracies 
of the levels have not been taken into account in 
making the figure. Only the group of levels close to 
the HOMO is shown, corresponding to MOs mainly 
due to the Mn 3d and Bi 6p states. Therefore, in the 
range of energies shown in Fig. [l^, there are twelve 
MOs for both spin orientations, and five additional 
MOs for the majority spin (only three levels are seen 
in Fig. [l^ due to the degeneracies) with main char- 
acter in the Mn 3c? states, which corroborates the 
d 5 configuration of the Mn atom also in the case of 
the free standing MnBi4 clusters. Only one of the 
Bi 6p MOs has a larger binding energy than the Mn 
3d states. This MO has a clear bonding configura- 
tion, with a p orbital from each Bi atom pointing 
in phase towards the central Mn, and includes some 
hybridization with the s symmetry orbitals of Mn. 
Above the energy levels shown in Fig. [l^ there is a 
'gap' of ~2 eV (~3 eV) for the minority (majority) 
spin, and well below the HOMO (~ —9 eV) four 
MOs coming from the Bi 6s states can be found. 

The most striking fact of the electronic structure 
presented in Fig. |l2| is the half metallic character it 
suggests: the minority spin valence states are all oc- 
cupied for both clusters, lying at least 0.2 eV below 
the HOMO. This situation holds for a large range of 
net charging of the MnBi4 clusters, showing the ro- 
bustness of the half metallicity. The most probable 
origin of this polarization of the Bi states is the Pauli 
exclusion principle: the majority spin MOs have to 
be orthogonal to the localized Mn 3d orbitals and 
the corresponding increase in kinetic energy pushes 
these levels to higher energies than those of the mi- 
nority spin states. As a result, the Bi 6p states are 
polarized oppositely to the Mn 3d moment. 

As a consequence of the half metallic electronic 
structure, the total spin of the tetrahedral MnBi 4 
clusters depends on the total number of valence elec- 
trons, as shown in Table [V. Taking into account the 
d 5 configuration obtained for the Mn atom, the MOs 
due to the 6p Bi states can still accomodate up to ten 
extra-electrons. In that case, with a net charge of 
-10 |e|, both majority and minority states are fully 
occupied and the magnetic moment of the cluster 
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TABLE IV. Total magnetic moment S (/Lts), as a 
functions of the net charge Q e (|e|) and total num- 
ber of valence electrons N e , for MnPbBi3, MnBi4 and, 
MnBisPo tetrahedral clusters with a central Mn atom. 
The cluster geometry was kept to that of the MnBi4 
tetrahedra in crystalline Ba^MnBin . 





MnPbBi?* - 


MnBi 


Qa- 
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MnBi 3 Po Qe " 


Qe 


N e 
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34 2 
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33 
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34 
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35 3 
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34 
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35 
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36 4 


-9 


35 
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36 
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37 5 


-10 


36 


4 


37 


5 




-11 


37 


5 









equals that of the Mn 3d shell (5/xb). If the net 
charge of the cluster is then reduced, the holes are 
always created in the majority spin MOs and, as a 
consequence, the total spin of the cluster is reduced. 
In particular, the formal charge of the MnBi 4 group 
in the manganites studied in this paper is -9 |ej (see 



section VII), which leaves one unocuppied majority- 
spin state. Therefore, consistent with the electronic 
structure of the clusters and the formal valence, we 
can expect a magnetic moment of ~4 /is on each 
MnBi^ - tetrahedron, as has been obtained in our 
calculations for the solids. The data in Table |y| also 
show that, if the charge transfer to the Mn tetrahe- 
dra is lower than 9 electrons, even a larger reduction 
of the magnetic moment could be expected. The ro- 
bustness of the dependence of the total spin on the 
number of electrons is also shown by substituting one 
of the Bi atoms by Pb, which provides one electron 
less, and by a Po atom, which provides one electron 
more than Bi to the cluster valence population. 

Since the rest of the atomic groups in the unit 
cell are formally closed shell, we can expect the 
main properties of the crystal to arise from MnBi4 
clusters, at least as far as the interaction— between 
MnBi^ - units is weak enough in the solidcS. With 
such simplified model we can already anticipate the 
following properties: 

i) the presence of a polarized hole localized in each 
MnBi 4 tetrahedron; 

ii) as a consequence of the reverse polarization of 
the valence band holes, the magnetic moment should 
be reduced considerably from the Mn 2+ value to 
near 4 \xb and, 

Hi) the FM order should exhibit a half-metallic 
band structure with approximately four unoccupied 
majority spin bands (one for each Mn atom in the 
unit cell); 

iv) the hopping of these holes between neighboring 
tetrahedra should be the main mechanism of electric 
transport. 



TABLE V. Population of the Mn d shell and mag- 
netic moments of different groups of atoms, as obtained 
from a Mulliken population analysis for FM Ba^MnBin , 
and the compounds obtained by substituting a Bi atom 
in one of the MnBi4 tetrahedra by a Pb and a Po 
atom, respectively. In all the cases the structure of the 
Bai4MnBin was maintained. The data correspond to 
the subtituted (MnBisX) groups, while those in brakets 
correspond the remaining MnBi4 tetrahedra. 





X=Pb 


Ba5 6 Mn 4 Bi43X 
X=Bi 


X=Po 


Q d 


5.14 (5.15) 


5.15 


5.16 (5.15) 


Md 


4.24 (4.33) 


4.32 


4.34 (4.32) 


l-lMn 


4.56 (4.63) 


4.62 


4.66 (4.62) 


HMnBi 3 X 


3.73 (4.39) 


4.35 


4.52 (4.32) 


f-^unitcell 


16.98 


17.59 


17.04 



These results are in fact corroborated by the cal- 
culations for the solids presented in the previous sec- 
tions, that take into account the full interconnected- 
ness of the lattice. 

In addition, we have performed calculations for 
some fictitious materials that provide further sup- 
port to the idea that an argumentation based on 
the MOs of the isolated clusters is still meaningful 
for the crystal. In one case we have substituted a 
Bi atom by Pb in one of the MnBi4 groups of the 
unit cell of Ba 14 MnBi n . According to Table IV the 
magnetic moment of the resulting MnPbBig - tetra- 
hedron should decrease to 3 from the 4 [Ib of 
the pure Bi cluster. In Tabic [v] we can see that, 
in the solid, this substitution is in fact accompanied 
by a decrease of both the total (by 0.61 \ib) and 
the local magnetic moment, which becomes 3.73 [Ib 
for the tetrahedron where the substitution has been 
performed while is almost unchanged for the rest 
of the tetrahedra (4.39 \lb)- The case of the Po 
substitution is a little bit more complicated. The 
magnetic moment of MnBiaPo 9- should increase to 
5 [is- This is again reflected in the results for the 
solid but in this case the observed change is more 
modest, 4.52 fis versus 4.32 /xg (see Table [v]). This 
different behavior is related to two different effects. 
In the first place, the Po substitution adds one elec- 
tron to the tetrahedron which is probably accom- 
panied by a derealization of the associated MOs 
(Wannier functions in the solid), thus the picture 
solely based on the electronic structure of the iso- 
lated clusters becomes less accurate. Another ef- 
fect comes from the fact that the band structure 
of Bai 4 MnBi 11; while very close to, is not totally 
half metallic. With the addition of one electron the 
half metallicity is reached, forcing the total magnetic 
moment to have integer values, which consequently 
diminishes it from 17.59 to 17 /is with the Po substi- 
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tution. As expected, all the changes of the magnetic 
moment observed upon chemical substitution occur 
without much modification of either the charge or 
the moment associated to the Mn atom, or its d 
shell. 



X. CONCLUSIONS 

In this paper we have presented and analyzed the 
results of first-principles, full-unit-cell, LSDA calcu- 
lations of CauMnBin and Ba^MnBin. The gen- 
eral agreement between our results and the available 
experimental information supports a microscopic 
framework within which to understand these com- 
plex materials, which show the simultaneous pres- 
ence of localized Mn magnetic moments and metal- 
licity. We also provide evidence for new properties, 
not yet corroborated by measurements, such as near 
half-metallicity of the FM manganites. 

We summarize our main conclusions: 

(i) Electron counting based on the Zintl pic- 
ture, which visualizes these materials as an ionically 
bound collection of dipositive alkaline earth ions 
and negatively charged ions and covalently bonded 
polyatomic anions (MnBi4 9 ~, B3 7 ~, and Bi 3_ ), is a 
good starting point for the description of the elec- 
tronic structure_of these compounds, as suggested 
by Gallup et alM 

(it) The Mn ions should be regarded as divalent 
Mn 2+ , rather than the original picture that consid- 
ers a d A configuration (3+ valence) for the Mn ion. 
Charge balance is maintained by the introduction of 
one hole into the Bi 6p states in the MnBi4 tetra- 
hedron; a heuristic connection Ga 3+ (Bi4) 12 ~ < ► 

Mn 2+ (Bi4) n ~ can be drawn, which also relates 
the charge states in the tetrahedron to the insu- 
lating or conducting character of the solid. The 
hole in the Bi4 tetrahedron can hop to neighboring 
MnBi4 units, accounting for the (hole) conduction 
in Cai4MnBin and Ba^MnBin. 

(Hi) In spite of the majority Mn 3d states being 
fully occupied, the experimental magnetic moment 
of ~4 /iB per formula unit (i.e. per Mn atom) is 
recovered in our calculations, since the hole in the 
Bi4 tetrahedron is parallel to the Mn moment (the 
unpaired electron is antiparallcl) . 

(iv) The FM compounds are close to half metal- 
lic, with the minority bands almost completely oc- 
cupied. 

(v) The AeMPn crystal structure introduces its 
own peculiar magnetic coupling. The MnBi4 tetra- 
hedra (the magnetic units of these compounds) are 
distributed in two disjoint interpenetrating three di- 
mensional networks. Magnetic coupling is ferromag- 
netic within each network, and AFM order arises if 



there are AFM interactions between tetrahedra from 
different networks. These anistropic magnetic cou- 
plings cannot be understood within a traditional, 
isotropic RKKY scheme. 

(vi) Our calculations also provide some insight 
into the puzzle of how the transition temperatures 
can be as high as 70 K for this family of compounds 
when the Mn-Mn are more than 10 A apart. One 
point is that the Mn-Mn distance is not the mean- 
ingful parameter, rather the relevant parameter is 
the Bi2-Bi2 distance which determines the hopping 
along the -Mn-Bi2-Bi2-Mn- chains. 

(vii) We speculate that full three dimensional or- 
dering is enhanced by the comparatively strong FM 
coupling within each network, analogously to the sit- 
uation in layered cuprates such as La2Cu04 which 
order magnetically above room temperature in spite 
of very weak interlayer exchange couplings. 

Our new picture of Mn in the AeMPn struc- 
ture should impact the interpretation of the va- 
lence fluctuation behavior of YbuZnSbn, where in- 
dications, are that both Yb 2+ and Yb 3+ ions are 
presenta. Dipositive Zn destroys the electron count 
that would result in semiconducting behavior (viz. 
Cai4Ga 3+ Asn). The system accommodates by 
drawing an electron from the Yb ions, which is pos- 
sible because of the near degeneracy of the Yb 2+ and 
Yb 3+ configurations. However, it may also be possi- 
ble to draw an electron from the 5p cluster orbitals 
of some of the Sb4 tetrahedra, as is done in every Bi4 
tetrahedron in Cai4MnBin and Ba^MnBin, with 
an accompanying unbalanced spin. Further study 
will be required to develop a consistent picture of 
this fluctuating valence system. 

Our description also gives new insight into the in- 
terplay between the conduction and magnetic prop- 
erties, and may play a crucial role in the CMR be- 
havior observed in some of the 14-1-11 manganites. 
In particular, it provides a simple picture of the pres- 
sure induced-, ferromagnetic transition observed in 
Sri4MnAsnEJ, which seems to be accompanied by a 
semicondutor-metal transition, and was previously 
interpreted as signature of a RKKY type of mag- 
netic interaction. Compressing the AeMPn struc- 
ture causes an increase in the interation between 
the MnBi4 units, with the corresponding increase 
of the hopping probability of the localized valence 
band holes. These holes are simultaneously respon- 
sibles of the conduction process and vehicles of the 
exchange interaction between Mn atoms and, there- 
fore, a decrease in the resistivity of these systems is 
necessarily accompanied by stronger magnetic inter- 
actions. 
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